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I. INTRODUCTION 


1.1 Background 

Continuous monitoring of ocean surfaces, utilizing the remote 
sensing capability of microwave radiometers, may provide Important 
Information applicable In the fields of oceanography, meteorology, 
marine transportation, and other related areas. Ocean surface 
characteristics such as salinity, molecular temperature, roughness, 
and sea Ice content can be determined from electromagnetic measure- 
ments, as reviewed by Tomlyasu [1]. These properties are parameters 
of the function describing the brightness temperature of the water, 
which Is also dependent on the frequency, polarization, and viewing 
angle of the system. The radiometer however, does not measure the 
brightness temperature directly. Instead, It measures the antenna 
temperature, which Is the brightness temperature ''smoothed'* by the 
system weighting function (antenna pattern). The blurring of the 
brightness temperature distribution Is experienced by all realizable 
antennas. 

To calculate the water brightness temperature from the measured 
antenna temperature, an antenna desmoothing problem Is encountered. 
This Is a classical problem In the field of radio astronomy [2]-[4], 
and some of the earlier methods used to solve this problem will be 
adapted to the remote sensing radiometer system under Investigation. 
However, radio astronomers have been able In recent years to reduce 
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the smoothing effect of the antenna by using highly efficient 
narrow beam aerials. This allows the observed antenna temperature 
distribution to closely approximate the brightness temperature pro- 
file* thus eliminating the desmoothing problem. To achieve this 
high resolution In the antenna pattern at microwave frequencies. It 
Is necessary to construct aerials with relatively large physical 
dimensions. However* the radiometer systems used for continuous 
ocean monitoring are alr-bome* placing considerable size limi- 
tations on the antenna. It may not even be desirable to use narrow 
beam aerials for fast moving air-borne radiometers* since half of 
the main beam must remain In the target area throughout the 
measuring period [5]. In any case* desmoothing methods must be 
applied to the wide beam antennas used In the remote sensing system. 

Desrooothlng of the antenna temperature distribution Is 
dependent on the successful calculation of a unique solution for a 
Fredholm Integral equation of the first kind. This Integral 
equation occurs In other Indirect measurement problems such as 
filter spectrometry* Inference of structure from diffraction 
patterns, etc. [6]-[8]. In every case the Fredholm equation Is 
highly unstable and yields unrealistic results when the measurements 
are known only with moderate accuracy (111-posed). Several matrix 
techniques have been developed [8]-[ll], along with the necessary 
smoothing procedures* which will Invert the Integral. For this 
antenna desmoothing problem* a more efficient Fourier transform 



Inversion method was developed and was then applied to an experi- 
mental finite wave tank system, taking Into account the temperature 
contributions from the surrounding sources (earth and sky). 

With the advent of the Fast Fourier Transform (FFT) algorithm 
[12],[13], the discrete frequency spectrums of various functions are 
comi'^uted very efficiently and economically. By placing the antenna 
temperature equation (Fredholm equation) Into a crosscorrelatlon 
for^ii, Fourier transform methods utilizing the FFT algorithm, allows 
the Inversion of the Integral to proceed efficiently. Also, since 
the frequency spectrum of the power pattern Is made available by the 
FFT, a better understanding of the smoothing effect of the i.er1a1 
may be obtained from examination of Its spectrum. 

1.2 Objective 

Microwave radiometer antennas tend to blur and smooth the 
brightness temperature profile of a target due to the coarse main 
beam and side-lobe level of the weighting function (power pattern). 
To recover the true brightness temperature distribution, FFT tech- 
niques are used to Invert the antenna temperature measurements. It 
Is the purpose of this research to 

a. verify the validity of the two-dimensional modeling, reported In 
reference [14], by comparing It with three-dimensional compu- 
tations. 

b. determine the stability and sensitivity of the Fourier Inversion 


method to profiles of errors In the antenne temperature. 

c. Isolate the causes of Instability and Illustrate them with 
examples. 

d. recommend and Implement modifications to smooth (filter) the 
Inherent Instabilities In the Inversion. 

e. Investigate the stability and sensitivity of the modified method. 

f. apply the modified techniques to actual measurements. 

II. THEORY 

2.1 Brightness Temperature Modeling 

All bodies that maintain a molecular temperature above 
absolute zero radiate electromagnetic energy that can be measured 
as a brightness temperature at a specific frequency. If the 
radiating object Is not a black body, then the brightness temperature 
will be less than that of an Ideal radiator. For opaque bodies, the 
brightness temperature can be expressed from Kirchoff's law [15] as 

Tb-eT„ 0) 

in which T^ Is the brightness temperature, e Is defined as the 
ability of a body to radiate energy (emissivlty), and T^, Is the 
molecular temperature. Since water Is considered to be opaque. Its 
brightness temperature will be given by (1). 

When (1) Is used to describe the brightness temperature of 
water, two polarizations of Its emissivlty must be considered. 

These polarized emissivltles can be described by the Fresnel 


equations for smooth flat water surfaces as 
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(2a) 


where 


e cos 0 


-/c - sin* 
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Th- 


COS 


0 - /e - s1n*0 


(2b) 


e cos 0 We - s1n*0 cos 0 + Jt - s1n*0 

and where 0 Is the Incidence angle of the energy, t Is the complex 
relative d1e1ec1;ic permittivity of the water, and v,h Indicate 
vertical and horizontal polarizations, respectively. Substituting 
(2b) Into (2a) for the vertical (horizontal) polarization, the 
corresponding emissivltles can be written as 


4 p COS0 




where 


4(pc' + qe")cos0 
(e‘ COS0 p)* + (e" COS0 + q)* (cos0 + p)* + q* 


p • {F cosy 
q ■ /F slny 

r » \/c(e' - s1n*0)* + (e”)*] 


(3a) 


(3b) 


Y » 5 tan’* (~ 


) 


2 e' - s1n‘0 

In which e' and c" are the real and Imaginary parts of the complex 
relative permittivity (c ■ c' - jE"). These equations are derived 
In Appendix I assuming that the air-water Interface Is a smooth flat 


surface. 
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To complete, the water's emlsslvlty model, the Debye equation, 
as modified by Stogryn [16], Is assumed for the complex permi- 
ttivities of (3). This equation Is 


I 

e 



e. + 


^0 • £«> 

1 - j2VTf 


+ j 


o 

airejf 


(4) 


where 

Eg • static dielectric constant of the solvent 

ca» * high frequency dielectric constant of the solvent (4.9) 

c* • permittivity of free space (farads/ni) 

a - Ionic conductivity of the dissolved salt (mhos/m) 

T ■ relaxation time (seconds) 

Stogryn has expressed the static dielectric constant e and the 
relaxation time t as semi -empirical equations which, based on a 
number of measurements, are written as 


6o (T,.M) • eo(Tm.O) a (N) (5) 

2rt(T„.N) - 2irc(T,.0) b (M,T«) (6) 

where 

eo(Tm.0)»87.74-0.40008 T„ + Q.aSSxlO-^Tm* + 1.410xl0“*T„* (7) 

a(N)»1.0-0.2551N + 5.151x10-* N* -6.889xlO-*N* (8) 

2itt(T||. 0)-1. 1109x10-*" -3.824x10-** Tm + 6.938xlO-**’Tm* 

- 5.096 x 10-*%* (9) 

b(N,T)-0.1463xl0-*IITm + 1.0-0.04896N - 0.02967N* (10) 

and 


T|q ■ water molecular temperature (0 < < 40^C) 

N * normality (0 < N 3) 

The normality In (5) and (6) can also be calculated from a seml- 
emplrlcal equation In terms of the salinity. This equation Is 

N«9[l. 707x10-* + 1.205x10** S ♦ 4.058x10“* S*] (11) 

where 

S « salinity In parts per thousand (0 < S < 260) 

The Ionic conductivity of sea water found In (4) can be formulated 
as a function of temperature and salinity. This equation Is semi- 
empirical In nature and has the form of 

®sea water ^^m»S) ■ ®sea water (25,5) exp (-Ac) (12) 

where 

«sea water (25.S) ■ $[0.182521-1.46192x10**5 + 2.69324x10**5* 

-1. 28205x10“ V] (13) 

C - 2.033x10“* + 1.266 x10“'’a + 2.464xlO“*A* - $[1,849x10-* 

-2.551xl0“*A + 2.551X10-V] (14) 

A-25-T„ (15) 

which Is valid for 0 s 5 < 40. Equations (3)-(15) require only the 
salinity, molecular temperature, and frequency to calculate accurate 
profiles of the vertically and horizontally polarized emissivltles 
of sea water. 

To get a more accurate representation of the water brightness 
temperature, equation (1) Is modified to include the reflection of 
the sky*s brightness temperature from the water surface. The 


modification of (1) Is 


O'Ey) ^bs (16) 

h h h 

In which Tbw Is the water brightness temperature* and Tb$ Is tha sky 
brightness temperature. This equation was formulated assuming that 
there were no contributions from any other atmospheric sources. 

To complete the modeling of all the brightness temperatures 
considered In the system* profiles must be selected for the earth 
and sky. For the lack of more reliable distributions* It will be 
assumed that the radiation emitted from the earth and sky Is 
unpolarized and for the earth will be fixed at a constant value of 
300^K. The brightness temperature of the sky Is best described [17] 
by the equation 

Tbs * Vf 0 - ® ) 

where 

faff 1.12 Tair - 50»K 
To • » In (1 - 3/Teff). (19) 

In which Tjir “ 284®K In all sky profiles modeled throughout the rest 
of this report. This model takes Into account atmospheric absorption 
by the two main contributions* oxygen and water vapor* and Is a 
reasonable good approximation for horizontally stratlsfled 
atmosphere. 

Equations (17)-(19) present sufficient Information for (16) to 
yield a fairly accurate water brightness temperature distribution. 

In Fig. 1 typical profiles of water brightness temperatures for 


(17) 

(18) 
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vertical and horizontal polarizations are plotted using the 
enissivlty modeling described by (2)-(15). This model was calcu* 
lated for the parameters Indicated In the figure and will be used 
throughout the Investigation. 

2.2 Description of Laboratory System 

Formulation of the antenna temperature equations describing 
the performance of the laboratory microwave radiometer system, 
shown pictorlally In Fig. 2 and geometrically In Fig. 3, Is simpli- 
fied by assuming a two-dimensional modeling. The system Is now 
available at the Flight Instrumentation Division of NASA, Langley 
Research Center, Hampton, Virginia. The square wave tank Is 14 feet 
wide and has the ability of controlling the molecular temperature, 
salinity, foam, and wave height of the water. The antenna and 
radiometer, located at x', y', z‘, of Fig. 3, are positioned above 
the wave tank of width w, and can travel along an arc of constant 
radius p. The radial distance from the displaced antenna system to 
the center of the tank forms an Incidence angle with a normal to the 
tank surface, referred to as the system rotation angle a. For each 
rotation angle a, the system has the mobility to perform a complete 
360° scan. The antenna's maximum radiation Is directed along the 
observation scan angle 3. The 6 variations are assumed to be 
negligible, thus forming a two-dimensional system the validity of 
which remains to be proven. 

For the particular experimental Installation under Investigation, 
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the total antenna temperature Is due to three different sources; 
namely the water, earth, and sky. Assuming no contributions from 
any other atmospheric sources, the total antenna temperature is then 
expressed as 


Ta (e) • Taw (3) ♦ (3) + (3) (20) 

where Ta is the total antenna temperature. Taw is the antenna 
temperature due to the water, Tae is the antenna temperature due to 
the earth, and Tas is the antenna temperature due to the sky. All 
of the inversion methods yet to be developed will be applied to 
Taw. which can be isolated easily from (20). 

2.3 Antenna Temperature Equations 

% 

For two-dimensional considerations, Fisher [14] has shown that 
each contribution of (20) can be written as 


Taw (O) ■ 


Tbw (0* + 0 + o) G (0‘) d 0' 


( 21 ) 


-01 
r-0i 


Tae (0) ■ 


I. . 


Tbe (0* + 0 + ®) 6 (0 ) <* 0 


fOi' 


Tbe (0* + 0 + ») ® ^22) 


02 
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T»s (8) • f Tbs (»■ + 8 + o) G (»•) d »' 


-Tr 

+ Tbs o' + B ♦ o) 6 o') d 0* (23) 

o“ 

,Zir 

(24) 


.2ir 

O') • P 0 ')/| P o') d 0 ‘ 


and the limits of Integration are found from the geometry of Fig. 3 
as 


0t»i 


I ,w/2 ♦ 0 4<na, 

*•" t p cos « J T “ 


(25) 





Where P(0') Is the power pattern of the aerial. G (0*) Is the 
normalized gain function of the antenna, TbM> Tbe* and Tb$ are the 
brightness temperatures of the water, earth, and sky, respectively. 
The brightness temperatures In (21)*(23) were described In section 
2.1. Each term In (21)>(23) was Intentionally placed In the form 
of a crosscorrelatlon and can be evaluated. In the Fourier trans- 
form domain, by a multiplication. From the crosscorrelatlon 
theorem, the antenna temperature equations (21)-(23) can be written 
In general as 


( 37 ) 


T,(u) • Tb(«) S*(»i) 

where Is the Fourier transform of Ta(B) defined by 

<o 

Ta(“) ■ I T,(S) e-J“8d8 (28) 

mm 

T|^(u) and G*(<d) are similarly the transform of T|)(0) and the con* 
Jugate of the transform of G(B)* respectively. The computational 
advantages of (27} are realized from the use of the FFT algorithm to 
calculate the direct and Inverse transforms of the functions Involved. 

A better understanding of the smoothing effect of the gain 
function Is obtained by examining Its spectrum's Influence on (27). 
This equation emphasizes the Importance of zeros In G*(<d). It Is 
apparent from (27) that for those frequencies for which 

G*((Dij) * 0 (29) 

the Information contained In T()((d|() Is lost In the smoothing 
process. This means that Ta((o) has no knowledge of T(,((i>) at 
no matter what the magnitude of T|)((i)|(). Even If G"^((d) only 
approaches zero, spectral Information Is severely attenuated. It Is 
reasonable to expect to be small In magnitude at some large 
u since continous functions, such as antenna patterns, have high 
frequency components that approach zero. 

Fig. 4 demonstrates the aerial smoothing effect In the trans- 
form domain for arbitrary functions. The point labeled may be 
described as being a cutoff frequency beyond which G(b)) Is at or near 
zero. Starting with T|,(u>) (Fig. 4a) and multiplying by G(w) 






(Fig. 4b) results In the spectrum of (Fig. 4c). This pro- 
cedure may be explained as being a filtering process In which the 
Input signal Tb(w) Is filtered by the system'Hi transfer function 
G((it), resulting In the ouput signal Ta(u)). The figure shows that 
the frequencies above u^. are completely lost and below reduced 
In magnitude according to the gain function spectrum variation. 

The end result of the two areas of spectral attenuation Is the 
smoothing of the brightness temperature distribution In the spatial 
domain. 

As was stated In the introduction, radio astrononers have been 
able to by-pass the antenna smoothing problem. They accomplish this 
by using highly efficient narrow beam aerials which allow the antenna 
temperature to approximate the brightness temperature. This approxl- 
rr:%,t1on Is accurate because the gain function of the high efficiency 
antenna approaches the Dirac delta function. If G(()) Is a delta 
function, Ta(u)) equals Tb(u) In (27). In this case, there Is no 
loss of spectral Information which results In the measurement of a 
highly resolved profile. However, the antennas used for remote 
sensing must meet rigid constraints that force the use of aerials 
with pattern frequency spectrums that decrease and approach zero at 
high frequencies. 

The existence of u^, that was discussed previously, places 
G((d) In the role of a low-pass filter that rejects the high 
frequency signals subjected to It. This Implies the possible 
existence of brightness temperature profiles that may be virtually 




undetectable, such as 

Tb(») - I expUi^9) (30) 

where the are arbitrary and are larger than a>g of Fig. 4. 

These temperature distributions may be called Invisible distri- 
butions, since they are not accepted and not recorded by the system. 
Although the presence of what (30) represents in nature is 
impossible, since it would require the measurement of negative 
absolute temperatures, it may however occur on top of a realizable 
detectable distribution such as T (0) forming 

Tb(P) - Tt0) + jj A;„ exp(j«„0) (31) 

The plausibility of (31) in water brightness temperature investi- 
gations is questionable, although its usage in the study of the 
stability of inversion methods is advantageous. 

2.4 Inversion of the Integral Equation 

Experiments have been devised to investigate the state of water 
by remote sensing techniques. For this investigation, the actual 
measurements would represent the antenna temperature, as given by 
(23), where the parameter of importance is the brightness 
temperature of the water which is part of the integrand of 
(21). It would then be important to be able to determine the 
functional variation within the Integral from a knowledge of the 
integral Itself (measurements). This process is known as inversion. 

Since (21) can be cast into a multiplication form in the 


transform domain, as shown by (27), T^(0) can then be determined by 


Tbw(B) ■ (32) 

where F~^ Is the Inverse Fourier transform and Taw(u) Is the 
antenna temperature due to the water In the frequency domain. It 
should be pointed out that Is obtained from (20) after the 
contributions from the earth (Tje) and sky (Tas)» calculated using 
(22) and (23) and a knowledge of and Tbs* have been subtracted 
from the total measured antenna temperature. The advantage of 
using (32) for Inverting (21) Is that the for all observation 
angles (6's) and one system position (a) are processed In one 
efficient calculation with the aid of the FFT algorithm. This 
makes the method very attractive In the evaluation of large 
quantities of data. 

2.5 Stability Investigation of the Inversion 

Calculating a Tg^ from a known Tbw by (27)and then Inverting 
this using (32) gives a Tb^ different from the original by only 
t 0.02^K, as was demonstrated In [14]. This would Indicate that 
(32) Is an accurate Inversion method. However, this method was 
examined only with exact T^i,. Therefore, It Is Important to 
Investigate the sensitivity and stability of (32) when TgM(3) Is 
known with only moderate accuracy, that Is to say when (21) 1s 
Ill-posed. 
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The Integral equation of (21) Is In a crosscorrelation form 
and Is a special case of the Fredholm Integral equation of the 
first kind. However, the stability of the more stringent general 
form of the Fredholm equation will be Inspected, and the con- 
clusions drawn from this Investigation will apply to the special 
case of (21). In the course of the stability Investigation, the 
following relationship Is required. 

11m 

R ♦ • 
a 

The proof of (33) Is given In [18] and Is repeated In Appendix II. 

The general form of the Fredholm Integral equation of the first 
kind Is 

.b 

g(y) « K(x,y) f(x) dx ; a<x<b (34) 

a 

where g(y) Is a known (measured) function, K(x,y) Is the kernel of 
the Integral equation, f(x) Is the unknown function, and a,b are 
constant limits of Integration. If for any A, Asln(Rx) Is added on 
to f(x) and R ® , then from (33), g(y) remains unchanged and (34) 
may be written as 

g(y) ■ I K(x,y) [f(x) + As1n(»x)] dx ; a ^ x $ b (35) 


K(x,y) s1n(Rx)dx - 0 ; a s x i b (33) 



However, if R if given sone finite large value rather than infinity, 
then an arbitrary small variable 6(y) Is added to g(y) forming 




g(y) + «(y) 


K(x,y) [AsInRx * f(x)] dx ; a s x s b 


(36) 


a 

If the left side of (36) represents measured data, that differs 
from the correct measurements g(y) by some small experimental error 
6(y|, the Inversion would then yield large amplitude, high frequency 
oscillations attached to the correct solution of f(x). This indi- 
cates that the solutions to the Fredholm Integral equation of the 
first kind are not unique when experimental errors are considered. 

The solutions of the Ill-posed Integral equation of the first 
kind will result In sporadic oscillations, regardless of the 
Inversion method used. This also applies to the Fourier Inversion 
technique of (32). The source of the Instability in the Fourier 
inversion becomes apparent If the spectrum of the gain function 
G(u) Is examined, as was done In section 3.2. 

To explain the Instability in the Fourier inversion, (32) Is 
rewritten as 


Tbw(B) - F-* 


- Taw(m) ^ Isliilin 
'■ 6*(w) G*(a>)'’ 


(37) 


where Te((d) Is the additive experimental error in the transform 
domain, and Taw'(«>) Is the spectrum of the error-free antenna 
temperature. For some <d near or above of Fig. 4b, the second 


term on the right side of (37) may be large because Te(u} Is 
amplified considerably from Its division by the small G(o>) compo- 
nents. In the extreme case where G(ii)|() Is exactly zero, the 
Inversion of the 111-posed (21) will blow up since T(jm( 3) will have 
components of Infinite magnitude for any trace of Te(u>|(). This 
means that the Inversion Inadvertently exaggerates the Am terms of 
the Invisible portion of (31). This Inversion method can not be used 
In practical applications and a modification to the technique roust 
be developed. 

2.6 Restoration Method 

The discussion on equation (37) showed that It Is Impossible to 
obtain an exact solution using (32) when T^i^ contains experimental 
errors. To reduce the oscillations of the solution, a smoothing 
procedure must be Included, as was done for the matrix Inversion 
techniques [88-[n]. This will require the filtering of the large 
high frequency error terms of (37). 

Writing 1/G*(u>) of (32) as l/l-Cl-G*(oj)], and then performing 
a series expansion [2] results In 

Tb«(6) • F-‘ a„M [1+(1-G*M) + (1-G*(»))' + (1-G*(«))* 

+....+ (1-G*(u))" +....]) (38) 

The Infinite series expansion of 1/G*(u>) converges provided that 
|l-G*(w)j < 1. For most antenna systems used In radlometry, their 
gain patterns are symmetrical, smooth varying functions which 


Insure that G*(u)) Is always real and In most cases positive. For 
some cases where G*(a>) Is negative, the amplitudes are very small 
compared to the positive values. Since G(u) has been normalized so 
that Its maximum value Is unity, convergence of 1/G'^(u) Insures the 
recovery of a unique Tbw only If Taw does not contain any error. 

The Inclusion of error In Taw leads to nonunique and oscillating 
solutions for Tbw since the convergence of 1/G*(u>) equates (38) to 

(37) . The oscillating forms of Tbw *»*e mainly attributed to the 
Inclusion of higher order terms in the series expansion of 1/G*(u>). 
Thus, Tbif(B) can be restored to the extent depending on the error 
present In Taw(<»). The oscillations can be minimized by properly 
truncating the series expansion, which In effect filters the high 
frequency error terms. In practice, only the first few terms of 

(38) are necessary, as will be shown later. 

The rewriting of (32) In a series form as in (38), places the 
equation In a convenient form to develop a smoothing procedure In 
the recovery of Tbw(B). The technique Involves a truncation of the 
Infinite series and then a multiplication of the Tj^v^(o)) by the 
truncated series. The inclusion of one, two, three, or more terms 
of the series can be Interpreted as letting the values of Tbvy 


be equal, respectively, to 


Tbwo “ Taw 

(39a) 

Tbwi • + (Tji^ - G*Tbwo) 

(39b) 


^bW2 * ^bwi * “ G*Tbwi) 


(39c) 


^bwn ■ ^bw (n-i) * Om * ®*^bw(n-»)^ 
where * Implies crosscorrelation. The altered Inversion procedure 

reduces to an Iterative method, as Indicated by (39a)>(39d),and will 
be referred to as restoration [2]. The second term In (39b)-(39d) Is 
a correction factor which Is added to the values of the previous 
restored brightness temperature to obtain the newly created 
function. It should be emphasized that this Iterative method allows 
restoration of the brightness temperature of the water for all 
viewing angles within the limits of the wave tank. 

The restoration method approximates well the correct 
solution from the truncation of (38). This filtering process not 
only reduces the oscillatory part of the solution, but It also 
smooths the true T^^^ where high frequency terms are required to 
describe the rapidly varying parts of the curve. However, the Tpw 
curve Is a fairly smooth function In Itself. The restoration 
method will smooth only the discontinuous edges of the wave tank 
profile or those temperatures near the horizon, depending on the 
size of the wave tank. This property of the restoration process 
will be demonstrated later. 
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III. COMPUTATIONS 

3.1 Calculation of the Antenna Temperature 

Computations for the Investigation of the Fourier Inversion 
technique will be conducted using the 5^ Half Power Beam Width 
(HPBW) antenna of Fig. 5, whose beam efficiency to the first null, 
assuming rotational symmetry. Is 99.97%. This beam efficiency Is 
typical of antenna patterns used In state-of-the-art radiometers 
nowadays [19]. To use the FFT algorithm developed by Cooley and 
Tukey [12] to perform a discrete Fourier transform, requires that 
2*^ sampling points be known (N Is an Integer). All functions used 
throughout this Investigation are sampled 256 times In one 
revel outlon (360®) of Its observation angle. Two hundred fifty six 

I 

points allow 128 harmonics of the fundamental frequency (fo zir) 
of the function to be calculated. Ninety (90) harmonics of the 
normalized spectrum of the 5® HPBW antenna are shown In Fig. 6. 
Although an exact can not be located in this figure, the spectral 
characteristics of Fig. 6 and Fig. 4b are equivalent. 

To prove the validity of the two-dimensional modeling of the 
wave tank system, T^^ of (21) was calculated from Tj^^ of Fig. 1 
using different methods. A three-dimensional model [20], using 
numerical Integration for computations, is compared to the two- 
dimensional Fourier transform technique in Table I. These results 
were obtained for various system rotation angles (a's) while the 
radiometer was looking toward the center of the wave tank (B*0). 







TABLE 1 


a = 20” 


a = 40^ 


a = 60° 


a = 800 


ANTENNA TEMPERATURES OF WATER (Taw) 

FOR VERTICAL AND HORIZONTAL POLARIZATIONS 

WATER ANTENNA w = 14 feet WATER ANTENNA 

TEMPERATURE (Tgw. P = 13 feet TEMPERATURE (Taw. °K) 


VERTICAL POLARIZATION 



WATER 

BRIGHTNESS 

TEMPERATURE 

(TbW, OK) 


109.10 


114.44 


133.19 


176.88 


269.70 


3-D 

NUMERICAL 

INTEGRATION 

■ 

2-D 

FAST FOURIER 
TRANSFORMS 

109.09 

109.23 

114.48 

114.59 

133.40 

133.42 

176.86 

176.93 

231.32 

232.84 


HORIZONTAL POLARIZATION 
WATER 3-D 


(TbW. 


109.10 


104.16 


89.26 


64.71 


37.97 


3-D 2-0 

NUMERICAL FAST FOURIER 
INTEGRATION TRANSFORMS 


109.09 

108.96 

104.14 

104.01 

89.19 

89.07 

64.52 

64.51 

34.24 

34.17 

































As Talle I Indicates* the two methods agree favorably and the T^m's 
are nearly the same as the corresponding for Incidence angles 
up to about 60^. For radiometer systems with less efficient 
radiation characteristics and dominant crosspolarizatlon components, 
larger differences between brightness and antenna temperatures will 
be found for all Incidence angles and the two-dimensional modeling 
will not then be as valid. 

3.2 Sensitivity Investigation of the Fourier Inversion and 

Restoration Methods 

As was explained previously using (37), the Instability In the 
Fourier Inversion Is dependent on the spectrums of the error 
and gain function (G((i))). To demonstrate this, three sinusodlal 
errors of different frequencies (10th, 30th, and 50th harmonics) and 
of constant magnitude (1°K) were Individually examined. Each of the 
sine errors were added on to the error-free vertically and hori- 
zontally polarized water antenna temperature (T^^) calculated using 
the water brightness temperature of Fig. 1 . The spectrums of the 
errors Te(S) are spikes placed at their corresponding frequencies, 
as shown In Fig. 7. Even though the magnitude of T 0 ( 3 ) Is the same 
In each case, the error In the Inverted T|,|^ Increases as the 
frequency of the error term Increases (Figs. 8 and 9). This Is 
expected from examination of Fig. 6 which shows G(u>) decreasing as 
0 ) Increases, causing the higher frequency error to be amplified In 
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Fig. 9. Empirical and inverted (with sinusoidal 
errors) horizontally polarized brightness 
temperature profiles. 


the Inversion process* since It Is divided by smaller numbers. The 
amplification of high frequency error terms Is true for all wide 
beam antennas (narrow frequency spectrum), as was explained In 
section 4.2. 

Another type of error that can be examined Is that which occurs 
when only one out of the two hundred fifty six measurements Is 
Incorrect. This Is Impulse error and has spectral characteristics 
similar to white noise (Fig. 10). The spectrum of the Impulse error 
Is constant and extends to all frequencies. Therefore, the osci- 
llations that result from the division of this Ta(w) by G*(w) are 
sporadic In nature and are large In magnitude. This Is demonstrated 
In Figs. 11 and 12 whose curves were calculated from a T^^ that 
contained a 10^ Impulse error at the edge of the wave tank. Since 
these Inversions are completely unacceptable, the restoration 
technique was then applied to the containing the Impulse error 
(Fig. 13). This method yields an accurate T^v^ curve out to an 
Incidence angle of approximately l?o using only three restorations 
(four terms of (38)). The discrepancies, that occur at Incidence 
angles greater than 17®, are due mostly to the distortion of the 
true Tb^^ by the filtering effect of the restoration process. 

The sine and Impulse errors were Introduced to show the 
sensitivity of the Inversion method to the spectrum of the error, 
and are not representative of errors In actual measurements. An 
example of a more realistic error Is demonstrated by considering 
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Fig. 12. Empirical and inverted (with 10® impulse 
error) horizontally polarized brightness 
temperature profiles. 
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Fig. 13. Empirical and restored (with assumed 
injected impulse error) brightness 
temperature for vertical and horizontal 
polarization. 


the operating procedures for the measurements and the manner In 
which they will be utilized to perform the Inversion. To complete 
the Inversion, It Is necessary to have measured data for a complete 
revolution (360^). Since It was decided that a sampling of about 
every 1.4^ was sufficient, two hundred fifty six (256) points will 
be required for each Inversion. Making 256 measurements necessi- 
tates long periods of time, and It Is not a very practical pro - 
cedure. Measurements are usually made every 5° or 10^ and the 
others, necessary for the Inversion, can be inferred by interpolating 
between measured points. 

To see what errors are Introduced from the procedure described. 
In Fig. 14 the error-free Tj for vertical polarization was plotted 
every 5® (denoted by stars) for the first 40® and every 10® 
thereafter. A curve was then drawn through the points (*—*—*) to 
represent a profile of T^ for all observation angles. The same 
procedure was followed for the horizontal polarization with the 
results siiown in Fig. 15. The true Ta is also shown with these 
curves ( — -) and seems to indicate very little difference between 
the two. However, closer observation shows large vertical 
differences around the rapidly varying parts of the curves. These 
errors have considerable high frequency content In their spectrums, 
resulting In the unrealistic Inversions of Fig. 16 and 17. If part 
of the data (0®-40®) is from the solid (— -) curve of Fig. 14 and 
the remaining (40®- 180®) from the dashed (—*—♦—) curve of the 
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same figure, the Inversions are still undesirable as shown In 
Fig. 18. The same Is true for the horizontal polarization as shown 
In Fig. 19 as obtained from the data of Fig. 15. 

Since complete Inversion, using (31), of the data shown In 
Figs. 14-19 resulted In quite unrealistic solutions, the restoration 
procedure of (38) was used to smooth the Instabilities. Performing 
only three restorations on the data used to calculate Figs. 16-19, 
the restored vertical and horizontal polarizations of Tbw are 
compared to the true Tbw Figs. 20 and 21. When Ta Is measured 
every 5® over the first 40®, discrepancies In the restored Tbw 
occur at angles near the edge of the wave tank because of the d1ff1« 
culty In measuring the Ta around this rapidly varying part of the 
profile. If the correct Ta Is used over the wave tank (every 1.4® 
over the 0®-40® range), then less disagreement Is observed between 
the true and restored Tbw difference between the curves In 
the second case Is not due to the error In Ta completely. Instead, 
as was discussed previously, this discrepancy results from the 
filtering process of the restoration which attenuates the high 
frequency terms In the solution. 

All of the high frequency components, made available by the FFT 
algorithm, are required to describe the discontinuous edge of the 
wave tank's brightness temperature profile. By eliminating some of 
these terms In the restoration process, a smoothed edge is formed. 
The discontinuous edge can be restored If more terms In (38) are 
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Fig. 18. Empirical and inverted (using data from Fig. 14) water 
brightness temperatures for vertical polarization. 
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Fig. 20. Computed and restored brightness temperature 
profiles of water for vertical polarization. 
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Fig, 21. Computed and restored brightness temperature 
profiles of water for horizontal polarization. 


considerecl. but this would enhance the errors In and would 
approach the unstable Inversions of Figs. 16-19. It was observed 
that three restorations is usually the limit for most measured data, 
which was also concluded In [3] and [21]. 

3.3 Refined Measuring Procedure 

To calculate the water's brightness temperature profile for 
large Incidence angles using a finite size wave tank, the system's 
rotation angle (a) roust be varied. Since only a small number of 
restorations will normally be performed on measurements, the usable 
portion of each restored profile, for a given antenna position 
along the arc, will be limited. Therefore, to maximize the benefits 
of the restoration process, the selection of the best first esti- 
mate for Ti)m of (39a) must be examined. 

Performing three restorations on measurements permits only 
three corrections to the first estimate of Tbw in (39). This raises 
the question of whether T^ might be a better estimate of T^y^ rather 
than Tqm, in which case its use would result in improved solutions 
from the limited number of restorations performed. To aid In the 
selection of the term to be placed In (39a), the curves in Figs. 22, 
23 (p ■ 13 feet) and Figs. 24,25 (p ■ 26 feet), were calculated by 
restoring the error-free Taw These plots can be compared to the 
curves of Figs. 26-29 which were calculated by restoring the error- 
free Ta. These figures Indicate that for small a the accuracy and 
range of the restored profiles are approximately the same whether 
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. 23. Restored brightness temperature profiles of Taw ^or odd o and p = 13 feet. 
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Fig. 25. Restored brightness temperature profiles of Tg^ for odd a and p = 26 feet. 
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Fiq. 27. Restored brightness temperature profiles of Tg for odd o and p = 13 feet. 
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Fig. 29. Restored brightness temperature profiles of T- for odd a and p = 26 feet. 
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the restoration Is being performed on the T,'s or T^w's. However, 
idien large a are considered. It can be concluded that the resto- 
ration should be performed on the Ta’s for vertical polarizations 
and Tail's for horizontal polarizations. This occurs because both 
TbWy Ta curve up for large Incidence anglis, as Tbwj^ «nd Taw 
both curve down for the same angles which results In a better first 
estimate for T|^ of (39a). 

The profiles of Figs. 22-29 show the overlapping of each a 
curve, and the limits of the usable data may be set In each case. 
These profiles however were calculated from data taken at 1.4® 
Intervals. Since It Is not practical to make measurements at this 
Interval, a compromised measuring scheme was developed. It Is Indi- 
cated from [22] that measurements can be made every 2.8^ over the 
wave tank and every 5.6® over the rest of the observation angles. 
This routine would be repeated for all a from 0® to 80° In 10® 
steps. 

The measuring scheme described requires Interpolation to pro- 
duce the 256 temperatures necessary for the restoration technique. 
Straight line Interpolation will be used for all measurements. This 
type of Interpolation will Introduce more error Into the earth and 
sky measurements than It will for the wave tank measur«nents, since 
earth and sky data are taken at larger Intervals. However, Figs. 20 
and 21 show that the restoration technique Is not sensitive to 
errors over the earth and sky when only three restorations are 
performed. 


A second source of error must be considered when T^w Is Iso- 
lated In ( 20 ) for the restoration of horizontally polarized data. 

In actual measurements and T|>$ are not known. To obtain a good 
approximation to these profiles* the portion of over the earth 
and sky Is assumed for T |)0 and Tb$> respectively. This allows the 
evaluation of Tae *nd Tas used In (20). To Improve the approximate 
Tbe profile* It will be assumed that 

Tbe(0) rTa(-0i-12O) (40) 

for - 01 - 12 ® < 0 < *01 and 0i < 0 < 02 + 12®. It Is more reasonable 
to assume that Tbe remains constant value In this range 

rather than decreasing as Ta does* since the earth Is homogenous In 
the area near the wave tank's edge. 

The final two-d1mens1ona' restoration program, utilizing the 
measuring scheme discussed* Is shown as a subroutine In Appendix III 
If the appropriate Ta measurements are made using the Tbw's of Fig. 

1 and placed In the computer program, the Tbw profiles of Figs. 30 
and 31 can then be calculated. These curves show that for the 
system parameters considered, the restored Tbw profile Is accurate 
to within 1®K out to an Incidence angle of 67® for p « 26 feet. 
Although the calculated Tbw beyond these limits have less accuracy 
than 1 ®K* they are closer to the true Tb^ than any other data 
available. 

The large Incorrect TbM that occur In the calculations of 
Figs. 30 and 31 are rejected by the computer program, which accounts 
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and edited restored a profiles (****). 
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for the breaks In the restored profile. These discrepancies are 
due mainly to the Interpolation of the measured data, and are more 
prominent when a 26 foot antenna boom Is considered. This occurs 
because the range of 6 within the limits of the wave tank for the 
26 foot boom Is small, causing the rapidly varying antenna 
temperatures to be Interpolated with considerable error. The result 
Is a reduction In the usable portion of each a profile, producing 
the breaks In the edited Tbw curves shown. However, It Is possible 
to fit a curve through the points made available by the computer 
program, yielding a fairly accurate and continuous Tbw profile. 

The discrepancy In the restored Tpw for each a curve near the 
edge of the wave tank In Figs. 22-29 Is explained as being caused 
by the attenuation of the high frequency terms of the solution. 

When an Infinite wave tank Is considered, the high frequency terms 
are needed to describe the rapid changes In the Tbw at horizon. 
This Is shown In Fig. 32 where the restored Tbw 'is calculated from 
an error-free Taw* Since both vertical and horizontal polarizations 
of Tbw vary rapidly at the horizon, nothing is gained by restoring 
Ta Instead of Taw* Therefore, an accurate restored Tbw can be 
expected out to an Incidence angle of about 80^ for Infinite wave 
tank considerations, which Is representative of open sea 
measurements. 

The restoration method, demonstrated In this section, restores 
the water brightness ten^>erature with good computational efficiency. 
The computer time required to perform each additional restoration. 
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Fig. 32. Computed and restored brightness temperature profiles for an infinite wave 
tank geometry. 


once the Fourier transfonns of all pertinent functions have been 
found. Is on the order of 0.05 seconds. A total time of about 1-2 
seconds Is required to perform a complete restoration for one a and for 
one polarization. The computational advantages of this method are 
beneficial In the evaluation of large quantities of data. 

3.4 Restoration of Sea Water Measurements 

Radiometric measurements of sea water at 7.55 GHz were made at 
Cape Cod canal, as reported by Swift [23]. The E- and H-plane 
patterns of the antenna used for the measurements are given In Fig. 

33 with beam efficiencies of 80.34/69.34, 88.98/76.76, 95.68/88.94, 
97.06/92.66, 98.06/95.55, 99.81/99.83% for the E-/H- planes within 
the observation angles of 10®, 20®, 32®, 40®, 58®t and 108®, 
respectively [20]. The spectrum of the E- and H-plane patterns are . 
given In Figs. 34 and 35, and show some negative terms. These 
negative components are small In magnitude and will have little 
effect on the solution for the small number of restorations 
performed. 

Using Infinite wave tank geometry approximations, a T^m was 
cal r elated by performing three restorations on the measurements 
reported by Swift. The restored T|)m and the measured Ta ^re shown 
In Fig. 36. A very close agreement between the two Is Indicated 
for Incidence angles up to about 75®# as would be expected for 
efficient antennas with negligible crosspolarizatlon components. 
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Near the edges of the tank there are large disagreements between 
the two because the sky contributions are beginning to become 
significant and the crosspolarizatlon contributions become more 
evident. The measured Ta profile Is not smooth, and exhibits flat 
spots due to uncontrollable parameters such as periodic surface 
roughness, sunlight reflections, and the close proximity of the 
Irregular shoreline. 

To Interpret the measurements and to shed some insight In the 
validity of Stogryn's mathematical model [16], 1n Fig. 37 the 
restored Tbw was plotted and compared to that obtained using 
equations (2)-(16) for the parameters indicated, which were measured 
Independently at the site. It Is clear that the restored values 
have the same profile characteristics as the ones computed from the 
mathematical model. There are some distinct disagreements at large 
incidence angles, but It should be emphasized that the restoration 
method is not valid in the region near the edges of the wave tank 
as was shown In Fig. 32. The computed values were based on a smooth 
surface mathematical model where the restored temperatures were 
based on measurements of sea water which may have had some surface 
Irregularities during the course of the experiment. 

IV. CONCLUSION 

The validity of the two-dimensional modeling used In this 
Investigation was verified for efficient radiometric antennas that I 
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Fig. 37. Computed water brightness temperature from empirical equations and restored 
profile from measurements at Cape Cod canal. 


contain negligible crosspolarizatlon components. This type of 
antenna Is used in systems employing state-of-the-art radiometers. 
The validity of the modeling was demonstrated by comparing the 
calculated antenna temperature of a three-dimensional numerical 
integration method with that of the two-dimensional Fourier trans- 
form method. The results of the two methods agreed favorably for 
efficient antennas and any differences may be attributed to the 
neglected e variations and crosspolarizatlon contributions in the 
two-dimensional system. 

It was shown that the complete Inversion of the two-dimensional 
antenna temperature equation is unstable and yields unrealistic 
solutions when experimental errors are considered. This is true 
regardless of the Inversion method used, as was proven using the 
general form of the Fredholm Integral equation of the first kind. 

The Instabilities appear in the Fourier inversion because of the 
division of error terms by small gain pattern components in the 
transform domain. This was demonstrated by inverting antenna 
temperatures that contained sine and Impulse errors. These examples 
showed considerable amplification of the high frequency error terms 
caused by the spectrum of the antenna pattern. 

By performing a series expansion of the Inversion equation, the 
inherent Instabilities may be filtered by properly truncating the 
series. This ability of the restoration technique was demonstrated 
using sampled data of typical measurements which contained error, 
and indicated an insensitivity to errors in the earth and sky data. 


Although the restoration technique filters the true brightness 
temperature as well as the Instabilities, a large portion of the 
water brightness temperature profile Is obtainable by editing the 
restorations of measurements made at different system positions in 
a finite wave tank geometry. For open sea measurements 88% of the 
entire profile can be accurately calculated using the restoration 
technique. 

The solution of the integral equation along with the appro- 
priate smoothing technique can be applied to large amounts of data 
with ease. The two-dimensional restoration was used to restore sea 
water measurements made In an infinite wave tank geometry. These 
restorations compared favorably with values obtained from a semi- 
emplrlcal mathematical modeling of a smooth surface. 

The Fast Fourier Transform algorithm has enabled the solution 
of the two-dimensional antenna temperature integral, formulated as 
a crosscorrelation, to proceed with the utmost efficiency and speed. 
To obtain an accuracy In the restoration method comparable with the 
measuring accuracy of the state-of-the-art radiometers (less than 
O.l^K), the 6 variations and crosspolarization contributions must 
be considered. This will require a three-dimensional Inversion 
which Is now under consideration. The advantages of the Fourier 
restoration method will be felt greatly in this circumstance. 
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Appendix I. Vertically and Horizontally Polarized Water 
Emissivities for Smooth Flat Air-Water 
Interfaces 
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The emissivity of water that has a smooth flat surface is 
expressed as 

e=l-|r|' (A-1) 

where e is the emissivity and T is the electromaqnetic reflection 
coefficient. The air-water interface is assumed to be abrupt 
(Fig. 38), which allovys the reflection coefficient to be described 
by the Fresnel equations as 

PqCOS^I-PxCOSP)' 

P 5. ^ 

^ riQCOS?}+n,cos(!) 

n,cos0-nocos0' 

^h ■ niCos|lfTiiCos0' 

where 

rijjS impedance of the air (free space) 
pj* impedance of the water 
0 = incidence anqle of the energy 
0 '= refracted anqle of the energy 
and v,h indicate the vertical and horizontal polarizations, 
respectively. 

Examining the vertically polarized reflection coefficient 
first, (A-2a) can be written in the form of 



(A-2a) 

(A-2b) 
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where 

Pq= permeability of the air (free space) 

Eo= permittivity of the air (free space) 

Pj= permeability of the vi«?ater 
El* permittivity of the water 

Since the permeability of the water (pj) is equal to that of free 
space (po). Snell's Law of refraction, which expresses the 
refracted angle in terms of the incidence angle, is given by 

sin*0' = sin^0 (A-4) 

This allows (A-3) to be written as 



(A-5) 


i cos(J+/l /a- sin*0 


El 


The permittivity of the water is described by the complex relative 

1' 

e. 


Ej 

permittivity, e = T“ » placing (A-5) in the form of 




E COS0 -/E-sin^0 


(A-6) 


E COS0 +ye-sin^f) 

E is a complex quantity given by 

e=E - JE (A-7) 

where e', e" are the real and imaginary parts of the complex 
relative permittivity. Therefore, 


/c-sin^^ = /r" e" JV 


(A-8) 


( 



1n which 
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sin'jD)' + {z'f 

Y = -I 

2 e'-sin^0 

Placing (Ar8) and {P<-7) into (A-6) results in 


Tv = 


e'cos0 - p+j(q-e”cos0) 
e'cos0 + p-j(q+c”cos0) 


(A-9) 


where 

p = /F cosy 
q = /F siny 

Making the denominator of (A-9) real gives 

{e')^cos 0 - r + (e“)^cos^0 + j(2e qcos0-2pe'‘cos0) (A-10) 

Fy = ■ ;; 

(e'cos0 + p)* + (q + e'cos0)^ 

Taking the square of the magnitude of (A-10) results in 



(e')'‘cos 0 - 2r(e') cos*0 + 2(e‘ )^(e'')^cos'*0 + r^ 
[(e'cos0 + p)* + (q + e"cos0)^]^ 

-2r(E")^Cus^0 + (e”) cos jl +4(e' )^q^cos^0 
[(e'cos0 + p)^ + (q + e"cos0)^]^ 


4p (e ) cos 0 - 8e e qp cos 0 

+ — I II I 111 . ■ 

[(e'cos0 + p)^ + (q +c''cos0)^]^ 

The vertically polarized emissivity can be written 
(A-11) as 


(A-ll) 
from (A-1) and 


4(pe' 4 qe")cos0[(e' )^cos^0 4- + 2e'pcos0 

[(e'cos0 4- p)^ 4 (q + e”cos0)^]^ 
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(e'')^cos (8 + 2e" q cos0 + q^] 
[(e'cos0 + p)* + (q e"cos(J)^‘]^ 
which reduces to the final form 


4{pe' + qe")cos0 


Ev = 


(A-12) 


(A-13) 


(e'cosU + p)^ + (q+e”cos0)^ 

Next, the horizontally polarized reflection coefficient is 

examined, and (A-2b) can be written as 


^h = 


1^ COS 0 - COS 0 ' 


'^1 1^0 
— “ COS 0 + /— COS 0 


(A-14) 


1 V '^0 

Using (A-4) and the fact that allows (A-14) to be placed in 
the form of 



(A-15) 


Th 


COS 0 


/T-: 

- ve-si 


2 

sin 0 


(A-16) 


COS0 + ^ £-sin 0 

where e is given by (A^7). Using (A-8), equation (A-16) takas the 
form of 


’’h = 


COS 0 - p+jq 
COS# + p-jq 


(A-17) 
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where p and q are given in (A-9). Making the denominator of (A-17) 
real gives 


Th ‘ 


cos 0-r + j2qcos?) 


(cos0 + p)^ + q^ 

Taking the square of the magnitude of (A-18) results in 


(A-18) 




2 2 2 
2 (cos fi-r) + 4qcos 


,2 . _ 2^2 


(A-19) 


[(cos^+p) + q ] 

The horizontally polarized emissivity can be written from (A-1) 
and (A-19) as 


^h 


Apcos’p + 4rcos*0 + 4prcosfl + 4cos^?>(p^ - q*) 


(A-20) 


[{COS0 + p)=^ + q^]- 

which reduces to the final form 

4pcosfll 

£. = _i 


(A-21) 


(cosU + p)* + 
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Appendix II. Integration with High Frequency Integrands 
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If K(x,y) is a piecewise continuous function for a < x ^ b. 


then 




Tim 

R CX) 


K(x»y)sinRxdx=0 


(A-22) 


To prove (A-22)» it can be assumed without loss of generality that 
K(x,y) is continuous. The interval of integration can always be 
broken up into a finite member of subintervals and proved for each 
of the subintervals where the function K(x,y) is continuous. The 
y variable in K(x,y) will be dropped since it is a constant with 
respect to the integration variable x. 

To begin the proof a change of variables is performed. 
Replacing x with t + ^ » <^esults in 

rb 

'x)sinRxdx = - 


K(3 


h ^ 

R 


K(t + — )sinRtdt 
R 


(A-23) 


Doubli.ngC A-23) gives 
.b 

K(x)sinRxdx " 


K(t)sinRtdt - 


h 1 

^■R 


ir 


K(t + j^)sinRtdt (A-24) 


a- I 
R 


Equation (A-,?4) can be exapnded as 

.b 




K(x)sinRxdx = 


IT. 


[K(t + ]^)-K(t)]sinRtdt 
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K(t)sinRtdt - 


R 


K(t+ f)sinRtdt (A-25) 




If K(x) is continuous for a s x s b, then R can be chosen sufficiently 

larqe so that |K(t + ^ )-K(t)| < 5/(b-a) for all t in the interval. 

Also, R can be picked large enough so that — < 6/2M, where |K(t)|<M 

R 

in the interval. Rewriting (A-25) gives 


rb 


2 | 


K(x)sinPxdx| 


.b- 


[K{t + )-K(t)]sinRtdt| 



rb 


K(x)sinRxdx| < fb-a] + 2 + 2 


whi ch reduces to 


,b 

I K(x)sinRxdx| < 6 


(A-27) 


(A-28) 


a 
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Since 6 is arbitrary, it can be set equal to zero in (A-28) 
forcing R -► «. 
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Appendix III. Program to Restore Finite Wave 
Tank Measurements 
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0014 
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0032 

003 3 
0034 
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0036 
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0040 
.0041 
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i 0047 
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370 FlmHAl (2F10.4) 
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IF 1 AXH ) .GT.0.01 AX(1)=-AX(I) 

IF (CX( .U .GT.0.0} CX(i)=-CX(I) 

•AX( n=10.0>!‘*lAX( 1 , 

372 CXU)-10.0=»4(Ca( 1 )/|'^.0) 

DO 37A 1=1 iblNnl 
IF (POL. Eu. 1.0) bxsm = Axm 
379 IF (ROL.EC.2.0) GXS(I)=CXm 
Du 875 I=81N91 ,biN21 
IF (PDL.EU.l.O) LXm=Axm 

375 IF (POL. EQ. 2.0) DXm = CX(l) 
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Fxn )-DX(24BiN9-l42> 
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FX( i+Bli^2 )=GXS( I ) 
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DO 377 1=1 f BIN 

. GX$( I )=FX( I ) . 

IF (i.L£.BiN2) 0X(1 )=FX( Ub IN2 ) 

377 IF (I.GT.BIN2) DXl U = FX( I-B IN2 ) 
call four (DX.FGXStBIN) 

Du 376 1 = 1 , BIN 
376 FXn )=1.0 

CALL FOUR (FX,BX,BIM 
DO 379 1=1 , BIN 

379 BXU )=TT4iiX( i)4CC<\|JG(FGXS(l ) ) 

CALL REVFUR (BX.BIN) 

DO 900 1=1 , BIN 
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0144 

0145 

0146 


OUT 
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0155 

0156 

0157 

0158 

015 9 
OUO 
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016 2 
0163 
016 4 
016 5 
0166 
0167 


TBES(J+BIN4) = TA( UBIN4I 


453 TAt 1 )=aX( 1 ) 
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♦ ^^OUTiNc TU uaLCULATE Trit AKTfcNNA TERPEKATURE DUE TO 

♦ THE WATER EKuM THE TOTAL ANTENNA TEHPERATURE 

IF IPuL.tU.2.) GO TG 92 
GU fU 93 

92 X=Pi/2.0-ALP«-DEL 

1 = 1 - 

80 IF ( ( G1N4-U-1)»0LL,LT.( THUU.O^FAOn GO TO 81 
IrUl 

GO TO 80 

81 DU »2 J=I,blN4 
Tbb5( J)=TA( J) . 

TbcS ( J+b IN4j=TAlJ+blN4I 

IF t CblN4"J + l)4bfcL.LT.(Thl»-l2.0^RAO) } TBES( J+B IN4) = TA( Ul 
IF I i biN4-J+lUUEL.LT.THU T 8 ES ( J + B I N4 » *0 .0 
TBESI J+B1N2)=TAU+BIN2) 

IF ( JX'OEL.Ll .(TH2+12*0-i‘RADn TBES ( J + B IN2 )*:TA ( U-BIN4 ) 

IF I J^DEL.LT.TH2} T BES ( J+Bl N2 ) =0 .0 
IF iJ^DEL.Gl.XJ TfitSI J+B1N2»=TA( J40IN2I 

82 TbESI J-i-3*l3IN4)^-TA( j4-3*iiIN4} 

CALL l-UUK (TBESfFFXtblN) 

Do 701 1 = 1 , BIN 

701 F»-X( I ) = TT»CGNJG(FGX6( I ) UFFX n ) 

CALL hcVFUR (FFXfBlN) 

DU 702 1 = 1, BIN 

702 TAU )=TA( I)-BtAL(FFX(I) ) 

♦ kUUiIiMt TfJ CALCULATE THE WATER BRIGHTNESS TEMPERATURE * 

♦ BY Trie METHOD OF SUCCESSIVE kESTOPATICNS ♦ 

A44 Vi ♦♦♦=!«♦ 4 Vi 4 ♦♦^ ♦♦♦♦♦♦ Jr ♦*♦♦♦ 4 ♦♦♦♦♦♦♦♦♦*♦♦♦♦ 44 

v3 CALL FOUR (TA,FfX,bINJ 
DO 5o0 1=1 , BIN 

500 bX( I I=FFX( 1 ) 

Du yOI .4=1 ,3 
OiJ 501 I=i,bIN 

501 FFXl n = FFX{ I )+(1.0-lT=i^CuNJG(PGXS( in )+vM4BX( n 
LALL kcVFJh (f-FX,bIN) 

4V4 4-ii445;-444 4ii‘ 44 4444-v >)c444444 4 444 4 4 4 4 4 4 4 4 444 444444 44 * * *44 444 

♦ FukMS LONTiNULUS 1 bW CUkVE FKO^* 9 ALPHA CURVES ♦ 

V 4:^t 4 4;;i 4 j)0!c 444 4«4444-4444 4 44 444 4 vs4 4 4 444 4 4 444444444444 44444 444 

I:> (kmO.Lc. 13.1) lNU=IiMDEL ( NPUN) 

IF ( .1H0.LE.a 3,1) STAR=NLST{ Nf^UN) 
ii' 1 K HU • LE • 1 3 • I i S I P= NL F(\ { ivi P LN ) 

IF (kha.GT.13.n STAk=NHSTFiNPUN) 

IF (KnO,GT.13.n STP=nHFi\(NRLN) 

IF (.XHO.GT.13.1) lND=iNOeHlNRUN) 


ii' 1 K HU • LE • 1 3 • 1 1 S I P= NL F(\ { ivi P LN ) 

IF (kha.GT.13.n STAk=NHSTFiNPUN) 

IF (KnO,GT.13.n STP=NHFi\(NkLN) 

IF (.XHO.GT.13.1) lND=iNOeHlNRUN) 

J =1 

0(1 BOO I = S1 ARrSTP 

mN 0L( n=( Ii\lJ+STMk-l30 + J)4DtL«DEG-*-ALPCtNRUN) 
J = J+1 

THW( i )=KEAL( FFX( iND + i n 
NRU!>i = NRUN+‘1 


GG Tu 6 
oQu J =1 

Lu 801 1 = 1 ,S1 P 
IF IT BW (II • i\ L • O • 0 ) 
IF (TbW( iJ .Nc.O.U) 
BOi IF ( TbW( I ) .Nb.0.0) 
NFT=J-1 
kcTURN 
END 


•j Bw( J )=T8W( 1 I 
ANGL ( jUANGLH ) 
J=.l+l 
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CAT6 = 74203 


i.OCOl 


0002 

0003 

0004 

0005 

0006 

000 7 
ocoa 
000=) 
0010 
0011 
0012 
001 ’ 
OOlH 

001 5 

0016 
0017 
OOIB 

001 9 


; bUbKUUTlNf F0UR<FUNfCuMP,.>jUt6) 

Si *>!«*♦**♦« >)5 >l«-^ * ♦ l«i ★ # 5)1 >!< Ilf ^ ★ 

♦ CUHP= FFI UF THE- GIVRN KcAE i-UNCTlCNf FUN ♦ 

sc ijc 4i ](: 4 St iti Hi ^ 4‘ 4i 4i « « 4 ai> i# »)i * 4 * 4 4 * 4i * 4( * 4> 4i 4i 4i ♦ « >!<; « « )i|( 

iNTcGER SLR 
CUi^MUN SCR.L 

COMPLEX £,CCMP<45bi ,LC(2bbl 

DIMENSION FUN(256I ,C(U9), REALL <5 1 2 1 , SCR (256 1 
fcWUi VALENCE (CC(1 ) iREALLili ) 

NUM2 = 2=<‘UUM 

ni 1 1 K = > . Nl IM^ . > 


C WU i V INL.C 

NUM2 = 2=<‘UUM 
Du 1 K=2tNUM2»2 
RLALl (K i=0 .0 

Du 2 K— lfi\UM2i2 
RcALLlK)=FUN(KJ) 

Lf I — / i A \ 


uu ^ 

RcALLlK)=FUN( rv (j # 

Ou i J=1 * NDM 
COMP ( JJ=^CC( J) 

CALL FrUUR(NUMf CCMP,NuM) 

UU A l~i»NUM 

COMi- ( i,) = CCMPU l/FLOATtNuM 
RETUki. 


:MP,NuM) 


4 : 4c sc St 4 : 4t Sr V ).'.v4t > 1 « S'St <= JftJit V* ScSc^t^tSt 414.' 414 ijcst ><=4 4 4 ♦>(ts)t>it5!t4' >#'***'!' 4*44444# =^>'4 

♦ LUMP- INVLkSi: FFT of THE GIVEN CORFLEX FUNCTlONf COMP 4 

*4444si44st SC44 4 v44444V444444 44 4 44*4 4 4X 44 444444 4444444 44**4*4 

ENTRY kEVFCh (CLMP,NUM) 

CmLL FrUUR tNUMfCUMPfNUM) 

NN=NUR/2 

CUMP( U=CGMP(NUM) 

DQ 5 L=1 ,NN 
t=CUMPI NUM-L+1) 

LUMP ( NUN- L + 1 )=CUMP'4L + 1) 

5 COMP( L+l)=c 

RETURN ■ 

END 
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CATE 
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*0 001 
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0 00.^ 
■001..4 
000“? 
OC06 
0 00 7 
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0009 

0010 

0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 
0 019 
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OOZl 
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0023 

0024 

0025 

00c 6 

^ 0 02 7 
0028 

0029 

0030 

0031 

0032 
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0034 

0035 

0036 

0037 

0038 

0039 

0040 
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004 2 
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0044 
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004d 
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1nTl.oc.K T , Ti i-iL ,THi a ,131U, oisR 
LQHi''uJN SLR,U 
LOMFLuX u,F,X( 25-.I 
DIMklMSiON SCK(21i.,) ,u3E(2 ) ,C( 129) 
BUUIVALEiVCE (£,C,S&(iJ) 

T1NC = '3IN/T 
LlNC=T/2 
L0C = LIIJC 
THTA=0 

cSb(iJ=C(THTA+l) 

CSc(^)=-L( 01N/4-i-l-’ 

LGC = LuC*-1 
LUC1 = L'JC“L INC 
E=E4(x<lOCII-XILCC) ) 

X( LUC1)=X( LOCI ) + X( LQt. ) 

X(LDC)=£ 

ThTA=TmTA+TINC 
IF {THTA-B1n/ 2)20,30, 3C 
lFtTHTA-6lN/4)8C,bO,5C 
CStl 1)=-CI falN/2+i-TriTA) 
CSE(2)=-C( THTA-IBIN/4-1 

» j : T ,'1 iL , \ 


■T'hTA) 


bU TO 60 

1 1 C / I 1 1 r 


T V ^-i / \ 


90 LUC=LbL+L INC 
GJ Tu 40 

91 1M2-LING) 92 ,93,94 

92 LliML = LlNC/2 
TINC=T1 nC+T 1NC 
Gu TO 10 

93 DO iuO LQC=2,T,2 
L L. o 1 = LOG "■ 1 
E=X(LUG1 )-X(LOC) 
X4 Lulx)=X(ICC1 )<-X(L0G) 

100 X(LijC)=£ 

94 CONI INUE 
DO. 9 J=1,T 
I — bvjR (J I 

IFl J~I)9,9,11 
11 F=XIJ) 

x( j)=x( n 

XU)|=F 
9 CJNTINUt 
HBTUf.N 
END 
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ABSTRACT 


Existing microwave radiometer technology now provides a 
suitable method for remote determination of the ocean surface's 
absolute brightness temperature. To extract the brightness 
temperature of the water from the antenna temperature equation, an 
unstable Fredholm integral equation of the first kind is solved. 
Fast Fourier Transform techniques are used to invert the integral 
after it is placed into a crosscorrelation form. Application and 
verification of the methods to a two-dimensional modeling of a 
laboratory wave tank system are included. The instability of the 
Fredholm equation is then demonstrated and a restoration procedure 
is included which smooths the resulting oscillations. V'ith the 
recent availability and advances of Fast Fourier Transform techni- 
ques, the method presented becomes very attractive in the evalu- 
ation of large quantities of data. Actual radiometric measurements 
of sea water are inverted using the restoration method, incor- 
porating the advantages of the Fast Fourier Transform algorithm for 
computations. 


